Influence of finite Hund rules and charge transfer on properties of Haldane systems. 
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We consider the Kondo-Hubbard model with ferromagnetic exchange coupling Jh , showing that it is 
an approximate effective model for late transition metal-0 linear systems. We study the dependence 
of the charge and spin gaps Ac, As, and several spin-spin correlation functions, including the hidden 
order parameter Z(tt), as functions of Jn/t and U/t, by numerical diagonalization of finite systems. 
Except for Z(n), all properties converge slowly to the strong-coupling limit. When Jn/t ~ 2 and 
U/t ~ 7 (the effective parameters that we obtain for Y2BaNi0s), As is roughly half of the value 
expected from a strong-coupling expansion. 



I. INTRODUCTION 

After Haldane's conjecture that integer-spin antifer- 
romagnetic Heisenberg chains should exhibit a gap in 
their excitation spectrunxl, there has been a,-caHsider- 
able amount of research in spin S — 1 systemsatJO. Fas- 
cinating aspects of these systems are the presence of free 
spin-lj/,2 excitations at the end of sufficiently long finite 
chainspL and the presence of a hidden string-topological 
orderoEl. 

The compound Y2BaNiOs is a candidate to a nearly 
ideal realization of the spin S — 1 antiferromagnetic 
Heisenberg chain and stimulated intense research on the 
system recently. The ID character is supported by ex- 
perimental evidence which shows that the exchange cou- 
plings transverse to the chains are vety small and unable 
to induce long range magnetic orderu. The representa- 
tion of the two S =1/2 holes per Ni by a S — 1 spin in 
the effective one-band model is, however, an intuitive but 
not clearly justified simplification, being the underlying 
assumption a large Hund rule acting on the two relevant 
Ni-orbitals. Such a point of view has been adopted in 
most theoretical works of the system, like Air alternative 
interpretation of specific-heat experimentsu in Zn-doped 
Y2BaNiC>5 which raised doubts about the presence o£ 
free spin-1/2 excitations near the end of long phainsEl. 
Theoretical work motivated by other experimentslilH for 
Ca-doped Y2BaNiOs, for which the NiO chains are not 
broken, but doped with holes also retaia p*jly the ground 
state of the 3d 8 configuration of Ni +2 llHrM 

On general physical grounds one expects the effective 
Ni-Ni hopping (via O) to be of the order^f the corre- 
sponding one in CuO chains (t ~ 0.85eVEJ). Further- 
more, spectroscopic data in atomic Ni show that the 
Hund rule leads to a ferromagnetic exchange J' H ~ 1.6eV. 
Thus, it seems that not only the triplet ground state of 
the 3d 8 configuration of Ni +2 , but also the excited singlet 
should be taken into account for a realistic description of 
the system. Except for a brief discussion an. the charge 
gap and Ni L3 x-ray absorption spectrun£3, this issue 
has been unexplored so far, to our knowledge. Even if 



the effective intratomic repulsion U is large compared 
with t, as we shall show, the properties of the system dif- 
fer from those of the strong-coupling limit. In particular, 
the expression A5 = 0.41049 J fotJiic spin gap in terms 
of the effective Ni-Ni exchange jtl'tlj, is no longer valid 
(at least if J is calculated perturbatively, see section IV). 
Thus, at least a qualitative study of the effects of a re- 
alistic Jh (instead of infinite) on the properties of these 
systems seems necessary. 

In this work we derive and study the Kondo-Hubbard 
model with ferromagnetic coupling, as an approximate 
(in the sense which will be clarified in the next section) ef- 
fective model for linear transition metal-oxygen systems, 
in which the relevant transition-metal orbitals are the e g 
ones. The model retains the effects of charge fluctuations 
at the transition- metal ions and finite Hund rules. Ex- 
plicit effective parameters are calculated for Y2BaNiC>5. 
Its version with antiferromagnetic coupling Jh has been 
extenshply studied in the context of heavy fermion 
systemsEil while the ferromagnetic case also adquired rel- 
evance in connection with the physics of ±he perovskite 
Mn oxides with giant magnetoresistanceE3, h eij j igj also 
closely related to the double-exchange model p?il and 
other models used to study these compoundsEfc . 

In section II, we explain the derivation of of the Kondo- 
Hubbard model as a low-energy effective model, in which 
some terms of lower magnitude, were neglected for sim- 
plicity. Section III contains the results for charge and 
spin gaps, spin expectation values and several spin-spin 
correlation functions, obtained by numerical diagonaliza- 
tion. Section IV contains the conclusions. 



II. THE LOW-ENERGY REDUCTION 
PROCEDURE 

In the simplest and most usual perovskite structures, 
the transition metal atoms are in sites of nearly cubic 
(Oh) symmetry, surrounded by six O atoms, lying in 
the directions ±x, ±y, ±z. In the particular case of 
Y2BaNiC>5, these octahedra are linked by their vertices 



and form well separated chains, making a nearly ideal 
one-dimensional compound. Since by far, the largest can.- 
tribution to crystal-field splitting is due to covalencytj, 
near the end of the 3d series, the few holes present in the 
3d shell of the transition metal enter the e g (d 3z 2_ r 2 or 
d x 2_ y 2) orbitals. This is due to their larger hybridization 
with the 2p(j orbitals (those pointing towards the transi- 
tion metal atom) of the nearest-neighbor O atoms. The 
starting multiband Hamiltonian for the system should in- 
clude the above mentioned orbitals and can be divided 
as follows: 

H = Hd + Hp + H p d + H pp , (1) 

where Hd (H p ) describes the on-site correlations inside 
the 3d (2p) shell. H p d contains the transition metal-0 
hoppings and repulsions. The 0-0 hopping is described 
by Hpp. The last three terms of Eq. ([j]) are a trivial ex- 
tension of similar terms extensively described and stud- 
ied in multiband models for the cupratesE2l and will be 
not reproduced here. The largest energies in the problem 
and the ingredients of the new physics when more than 
one relevant orbital per site is present, are contained in 
Hd- The Coulomb and exchange integrals among the 
e g spin-orbitals can be parameterized in terms of three 
Slater parameters Fq, F% and -F4, using usual methods in 
atomic physicsEl Denoting a[ a (b\ a ) the creation opera- 
tor for the aig d 3z 2_ r 2 (6 lg d x 2_ y 2) orbital at site i with 
spin cr, the result can be written in the form: 

H d = Ud^2{a^a t1 al l a ll + &J T &iT & !j, 6 a) 

i 

HUd-J' H )J2 a l a ™ b l' b ™' 

iaa' 
i a a' 

+(at T at i 6 ii 6 iT +if.c.)], (2) 

where U d = F + 4F 2 + 36F 4 and J' H = 8F 2 + 30F 4 . Be- 
cause of the neglect of the t 2g {d xy , d yz , d zx ) orbitals, 
Hd lost the invariance under rotations of the atom, but 
retains cubic (Oh) symmetry. From the two lowest exci- 
tation energies of atomic Ni, we obtain F 2 — 0.1600eV 
and F 4 = 0.0108eV, leading to J' H = 1.60eV. According 
to the theoretical interpretation of optical experiments 
in NiO, U d + J' H /2 = lOeVB The variation of J' H along 
the 3d series is only a few per cent, while the value of Ud 
is very similar to that calculated in the cuprates using 
constrained density-functional theoryEl More sensitive 
to the particular system are the hopping parameters (in- 
cluded in Hpp and H p d) and, particularly, the transition- 
metal to O charge transfer energy A (defined as the en- 
ergy necesary to take a hole from the ground state 3d 8 
configuration of the transition metal and put it in the 
p„ ocbital of the 2p shell of a nearest-neighbor oxygen 
atomEj), which increases to the left of the periodic table. 

To derive the effective Hamilto- 
nian for one-dimensional transition metal-0 systems, we 
employ the cell perturbation methodEfl. For simplicity in 
the explanation below, we choose the particular case of 
NiOe octahedra sharing O atoms along the z direction 



(present in Y2BaNiOs) and assume tetragonal symme- 
try. The p z orbitals of the O atoms lying between two 
Ni atoms are expressed irui.erms of Wannier functions 7r 
centered at each Ni atomtJ. Each NiOs cell, composed 
of the 3d orbitals at one Ni site, the 7r orbital at that 
site and those of the four nearest O atoms along the ±2 
and ±y directionscj, is solved exactly. To construct the 
low-energy effective Hamiltonian, only eight eigenstates 
of the cell are retained: for two holes, the Bi g triplet 
(which is essentially the ground-state of the Ni +2 config- 
uration plus corrections due to hybridization), and the 
first excited state, the B\ g singlet. For one and two holes 
in the cell, the lowest B\ g doublets are retained. These 
eight states are mapped into the corresponding ones of 
the Kondo-Hubbard model Hkh (Eq. |3|) at the corre- 
sponding site. The matrix elements of H in the restricted 
basis are calculated and mapped into the corresponding 
ones of Hkh- The effect of the remaining states of H 
could be included perturbatively but we neglect it for 
simplicity. To retain a simple and more general form 
of Hkh, we also neglect the dependence of the result- 
ing effective hopping on the occupation and spin of the 
sites involved. The resulting effective Hamiltonian has 
the form: 

Hkh = -t*^2 a ia a i+i<? + U a^a^alf a.^ 

i i 

— Jh ^2 S ia • Sib, (3) 

i 

where Si a , are the spin of the fermions represented 
by the hole creation operators a\ a and b ia respectively. 
These are effective operators with the same symmetry as 
those entering Eq. (||) , but which differ from them in the 
general case. In what follows a\ a and b\ a refer to these 
effective operators. 

The meaning of Hkh is easier to understand in the 
limit A> Ud S> J' H , tpd, where t p d is the Ni-0 hopping 
along the chain. In this case, the interactions U and J of 
the effective Hamiltonian (0) coincide with those of Hkh, 
i.e. U ~ Ud and Jh ~ J^j, while the effective hopping 
in ([}]) is given by the second-order process which carries 
a id^ z 2_ r 2 hole to the same orbital of a nearest-neighbor 
Ni atom: t = i 2 pd / 'A. However, the case A 3> Ud is not 
representative of charge transfer systems like Y2BaNi05, 
for which A < Ud- In this case, the states aLaJi &] CT |0) 
actually represent states with occupation close to one in 
the O 7r orbitals. 

As a consequence, U is mainly determined by A in- 
stead of Ud, and the hopping matrix elements become 
dependent on the occupation and spin of the two sites 
involvedEZl. As mentioned above, this dependence was ne- 
glected to keep a simple and more general form of Hkh- 

To estimate the parameters of the effective model for 
Y2BaNi05, we took the values of J' H and Ud mentioned 
above, and the (more uncertain) values of A and the dif- 
ferent hopping parameters in H p d and H pp were taken 
from work on NiOE3, with the p — d and p — p hopping 
parameters scaled with distance r as r~ 7 / 2 and r~ 2 re- 
spectively. The parameters of Hkh which result from 
the mapping procedure are U = 4.4eV, Jh — 1.2eV and 
t ~ 0.7eV. It is interesting to note that Jh has a very 



small sensitivity to the parameters of H. Instead, chang- 
ing A and the hoppings of H within reasonable values 
affects U by - 20% and t by - 30%. 

III. RESULTS. 

In this section, we study the behavior of the charge and 
spin gap, spin expectation values and spin-spin correla- 
tion functions of Hkh, using Lanczos diagonalizations 
in periodic rings of length L = 4, 6 and 8. The rapid 
increase of the Hilbert space with L prevents us to study 
longer even chains with the present state of the art, but 
as we shall show, some trends are already clear. The unit 
of energy is chosen as t = 1. 

A. Charge gap 

In Fig. 1 we represent the charge gap Ac = E(l) + 
E(—l) — 2E(0), where E(n) is the ground-state energy 
for n added holes to the stoichiometric system (which 
contains one a± g and one b\ g hole per site) . The result for 
L = 8 is compared to that of a polynomial extrapolation 
in l/L to estimate finite-size effects. These effects are 
small for U > 4. As expected, the gap increases with 
U and Jh- In the strong-coupling limit t = 0, the gap 
is Ao = U+ Jh/2. As t is turned on, but kept small, 
the leading correction to E(0) is of order t 2 /Ao, while 
those oi E(l), E(—l) are equal and of order t. Assuming 
a Neel background (alternating spin projections 1 and - 
1) the correction for one added or one removed hole can 
be calculated and is —y/2t. In both cases, it is more 
convenient to align ferromagnetically the spin at the site 
of the added or removed hole with those of its nearest 
neighbors. Thus, for large U , Jh, we estimate: 

A c = A - 2V2t. (4) 

In Fig. 2 we represent Ac — An as a function of U and 
Jh ■ The results agree with Eq. (||) in the strong-coupling 
limit. In the opposite limit, for U = 0, and small values 
of Jh, the results have important size effects, and a large 
positive value of Ac — A is not reasonable. However, 
the extrapolated results show a reasonable behavior and 
tend to small values in the limit of Jh = 0. In any case, 
the results for U > 4 seem reliable. From the parameter 
estimates for Y2BaNiOs given at the end of the previous 
section (U/t = 6.3, Jn/t — 1-7, t = 0.7eV), we obtain 
Ac ~ 3eV. This is sopewhat larger than the experi- 
mental value Ac ~ 2eVll3. This discrepancy is probably 
due to the fact that the charge transfer gap and possibly 
the hopping parameters cannot be transferred directly 
from Ref.Ej (which is a theoretical interpretation of op- 
tical spectra in NiO) to Y2BaNi05. 

B. Spin gap 

In Fig. 3, we show the spin gap As as a function of 
Jh for several values of U. For U = 0, the result for 
As has already been reportedE3. Here, for the smaller 



values of U and Jh, the finite-size effects are too large, 
and the extrapolated values to the thermodynamic limit 
(in some cases negative) are meaningless. However, for 
more realistic values of U, our results allow us to extract 
some conclusions. 

The qualitative behavior of As as a function of Jh was 
to be expected from the limiting cases: if Jh = 0, the 
model is equivalent to the Hubbard model plus L free 
spin-1/2 states, and therefore As = 0. For jLp — > 0, 
As ~ J H has been obtained using bosonizatiorc3. In the 
limit of large Jh, the low-energy physics of Hkh reduces 
to a spin-1 Heisenberg chain: 

t 2 t 2 

H Heis = jJ2Si-S i+1 , J= =WT7 -^, (5) 

i 

where Sj = S,; a + Sj&. This model has a spin gap 
Acp ls = c{L) J, where the constant c(L) depends on the 
size of the system. In particular c(4) = 1, c(6) = 0.72, 
c(8) = 0.59, c(oo) = 0.41049B The first three values 
of c(L) coincide within 1% with our results extrapolated 
to infinite Jh- As Jh increases, first As increases from 
zero, and as the system approaches the strong-coupling 
limit, As decreases with the effective spin-1 exchange J. 

In Fig. 4 we show As/Ag els as a function of Jh- For 
U > 4, the extrapolated values do not differ very much 
from those of L = 8. Note that even for large values of 
U, the ratio As/A-g els is considerably smaller than 1 if 
Jh ~ t. In particular, for the parameters estimated for 
Y 2 BaNi0 5 , A 5 /Af eis ~ 0.5. However, this ratio was 
assumed 1 to estimate the value of J from experimen- 
tal measurements of As- Using the same set of parame- 
ters, we obtain from the extrapolated valueSpAq ~ 240K, 
while the experimental value is As ~ 100KLrlL2l. In view 
of the approximations made in deriving Hkh, the un- 
certainties in the parameters of H, and the sensitivity 
of As and J to these parameters, the result is satisfac- 
tory. Part of the overestimate is due to ferromagnetic 
corrections to J in second-order in the intercell hop- 
ping, which involve virtual quadruplet three-hole states. 
These states are contained in H, but were projected 
out of the Hilbert space of Hkh- We have calculated 
J by the cell-perturbation method, including these cor- 
rections. The effective J is reduced from 0.098eV to 
0.088eV, but the result is quite sensitive to the param- 
eters of H. Copaparison with exact diagonalizations of 
a Ni20nclustei£j, shows that the second-order result of 
the cell-perturbation method is still an overestimation by 
a factor near 2, due to higher order corrections. 

C. Spin expectation values 

The behavior of the spin gap as a function of Jh and 
U displays a slow change from the weak to the strong- 
coupling regimes. In Fig. 5(a) we show the ratio of the 
spin gap to the effective exchange, as a function of U 
for L = 8 and different values of Jh- For any non- 
zero value of Jh, the strong-coupling limit is reached 
for sufficiently large U and the gap tends to the limit 
Ag ets = 0.59 J. Fig 5(b) shows the corresponding change 
in the total spin of both itinerant (&i g ) and localized 



holes as U is increased, in the lowest-energy state with 
total spin and projection S t = = I. In the limit of 
small Jh, both types of holes are decoupled and it is 
easier to flip a localized hole rather than an itinerant one 
from the S t = ground state. As a consequence, S bt = l, 
S* t = 1 — S bt = 0. In the opposite limit of very large Jh, 

the singlet states at each site ((aj+bjr — a lx^ii)l^)) can ^ e 
projected out of the relevant Hilbert space, and in this 
case, the following equality among spin operators at a 
given site can be proved: 2Si a = 2Sn, — Si. Summing 
over all sites: 2S at = 2Sfc t = S t . Thus, S% t changes from 
to 1/2 as Jh increases and S bt = 1 — S% t . The ef- 
fect of increasing U is to localize the itinerant a\ g holes 
and therefore, to contribute to the effect of Jh, reaching 
faster the strong-coupling limit. However, it is noticeable 
that the approach to this limit is very slow. Comparison 
of the quantities represented in Fig. 5 for different sizes 
suggests that this approach is even slower in the thermo- 
dynamic limit. 



states and charge fluctuations are allowed. The results 
of Fig. 7 are a hint that Z(ir) can be used as the corre- 
sponding parameter of a hidden order that also exists in 
the Kondo-Hubbard model Hkh- The difference in the 
behavior of Z(jr) as a funcion of U, Jh, in comparison 
with that of the other correlation functions, is an indi- 
cation that Z(tt) is more sensitive to the transition to 
the spin-gap state. That is precisely what one expects 
from a quantity playing the role of an order parame- 
ter. The other spin-spin correlations should be short- 
rangejJike. due to the opening of the spin gap for finite 
U, JffEJEl Bosonization results indicate that thiSjgap in- 
creases quadratically with the Hund coupling JhE3, and 
is thus very small for Jh — > 0. In this regime, the cor- 
relation length is larger than the maximum system size 
that we have studied and we are unable to identify the 
change in the behavior of the correlation functions. For 
larger values of Jh, the change of regime is captured by 
our results, and more clearly from the behavior of Z(w). 



D. Spin-spin correlations 

In addition to the spin-spin correlation functions 

#2(0 = ((S* a - S* b )(S* +la - S- +lb )), (6) 

we also study in this section the string correlation func- 
tion 

Z(j-i) = {S?eM™ £ W>- ( ? ) 
j=»+i 

This latter has been propposed as a hidden order param- 
eter for S — 1 chains to describe a hidden Z2X Z2 sym- 
metry breaking corresponding to the appearance of the 
Haldane gapQd. This symmetry has been first implicitly 
introduced ip an elegant variational approach for the ex- 
cited statesQ. In Fig. 6, we show the different correlation 
functions. 52(0 nas an on-site value 62(0) = 0.3 and for 
other distances 5^(0 < 0.02 for the parameters of Fig. 
6. The antiferromagnetic correlations are evident. They 
are larger for the localized holes than for the itinerant 
ones, as expected. For U = 0, a tenedency to antiferro- 
magnetic order with wave vector q — n is expected from 
the form cos(2fcp r) of the oscilating Rudcrman-Kittel- 
Kasuya-Yosida effective interaction between localized b\ g 
holes at a distance r mediated by the mobile a\ g ones 
with Fermi wave vector kp = tt/2. In the strong-coupling 
limit, the effective model Hn e is © leads to the same 
type of short-range correlations. 

The Fourier transform (Si(q) = J2i e ' iql , etc.) 
of some of these correlation functions for wave vector 
q = ir, is represented in Fig. 7. As U increases, Si(tt) 
and approach slowly the asymptotic value in the 

strong-coupling limit, as it was the case of the spin gap 
and spin expectation values already discussed. Instead, 
Z(jr) seems to saturate faster to a fixed value as Jh and 
U increase. 

In the strorus-coupling limit, Z(tt) is a signature of the 
Haldane stateatl. For our model, Z(tt) is a possible gen- 
eralization of this order parameter, when local singlet 



IV. SUMMARY AND DISCUSSION 

We have studied charge and spin gap, spin expecta- 
tion values, and several spin-spin correlation functions of 
a Kondo-Heisenberg model Hkh, for two particles per 
site. We have shown that the model can be considered 
as an approximate effective model for one-dimensional 
transition metal-0 systems, in which only the e g orbitals 
of the transition metals are relevant. Without any ad- 
justable parameters (taking the parameters of the orig- 
inal multiband Hamiltonian from NiO scaled appropri- 
ately with distance), we obtain from the effective Hkh, 
a charge and a spin gap of the correct order of magnitude 
for Y2BaNi05. The model has also been used as a sim- 
plified model for the manganites, and our results should 
be qualitatively valid in the limit in which all Mn ions 
are Mn+ 3 . 

For sufficiently large U or Jh, the charge gap of Hkh 
is approximately given by Ac = U + Jh/2 — 2v2£. The 
effective model Hkh contains the spin-1 antiferromag- 
netic Heisenberg model (also called Haldane chain) in 
the strong-coupling limit t <C Jh,U. We obtain however 
that for realistic parameters for Y 2 BaNi05 or when the 
effective Hund-rule exchange coupling Jh is not much 
larger than the effective hopping t, several properties dif- 
fer from the Haldane limit. In particular, the different 
dynamics of the itinerant and mobile e g holes (reflecting 
that they do not behave as part of the same spin-1 ob- 
ject) are clearly manifested in spin-spin correlation func- 
tions. In addition, the spin gap A5 is roughly half of 
that expected from a strong coupling expansions. This 
fact should be taken into account when the effective spin- 
1 exchange J is extracted from experimental values of A5 
and in the consistent interpretation of different thermo- 
dynamic experiments together with A5. One possible 
way to interpret our results for A5 when U ^> Jh ~ t, is 
that the effective spin-1 Heisenberg Hamiltonian Eq. (|^) 
is still valid, but higher order corrections in t reduce ap- 
preciably the second-order result J — t 2 / (U + Jh/2) for 
the effective exchange. Fourth-order corrections which 
include local singlet states as intermediate states are con- 



sistcnt with this reduction. However, when t ~ Jh, per- 
turbation theory ceases to be valid, and it seems more 
adequate to include the local singlets ((aL&t, — a\,b\y )|0) 
in our notation) in the model Hamiltonian. 
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FIGURE CAPTIONS 

Fig. 1. Charge gap as a function of Jh for several 
values of U indicated inside the figure. Solid symbols 
are the result for L = 8 sites. Open symbols correspond 
to extrapolation to the thermodynamic limit from the 
results of L — 4, 6 and 8, using a quadratic polynomial 
in 1/L. 



Fig. 2. Same as in Fig. 1 for the charge gap minus 
its strong-coupling value A = U + Jh/%- 

Fig. 3. Spin gap as a function of Jh for several 
values of U. Solid symbols correspond to L = 8 and 
open symbols to the extrapolated value. 

Fig. 4. Ratio of the spin gap As over its strong- 
coupling value (Af eis = c(L)J, where J = t 2 /(U + 
Jh/2), see text) as a function of Jh- The meaning of 
the symbols is the same as before. 

Fig. 5. (a) Ratio of the spin gap to the effective 
exchange Ag/J as a function of U for L — 8 and several 
values oi Jh- Jh = 2 (circles), Jh — 4 (squares), Jh = 6 
(diamonds) and J H = 20 (dashed line) . (b) z component 
of the total spin of the a\ g (bi g ) holes function 
of U for L = 8 and several values of Jh indicated by the 
same solid (open) symbols as above, in the lowest-energy 
state with total spin and z component St = S% = 1. 

Fig. 6. Spin-spin correlation functions (S£S? +l ) as 
a function of distance I for a\ g (S a (l)) holes, b\ g (Sb(l)) 
holes, the sum of both spins (S\(l)), and that defined by 
Eq. (1), for Jh = 2 and two values of U : U = 2 (circles), 
and U = 10 (diamonds). 

Fig. 7. Fourier transform of the correlation func- 
tions at momentum ir for the sum (Si(tt)) and difference 
(52(71")) of the spin of both types of holes at a given site, 
and the hidden order parameter Z(ir). 

Fig. 8. Size dependence of Si(q) and Z(q) for Jh = 2 
, U = 10. 
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